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Abstract 

In this work a series of methods are developed for understanding the Friedmann 
equation when it is beyond the reach of the Chebyshev theorem. First it will be 
demonstrated that every solution of the Friedmann equation admits a representation 
as a roulette such that information on the latter may be used to obtain that for the 
former. Next the Friedmann equation is integrated for a quadratic equation of state 
and for the Randall-Sundrum II universe, leading to a harvest of a rich collection 
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of new interesting phenomena. Finally an analytic method is used to isolate the 
asymptotic behavior of the solutions of the Friedmann equation, when the equation 
of state is of an extended form which renders the integration impossible, and to 
establish a universal exponential growth law. 

Keywords: Astrophysical fluid dynamics, cosmology with extra dimensions, 
alternatives to inflation, initial conditions and eternal universe, cosmological appli¬ 
cations of theories with extra dimensions, string theory and cosmology. 
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1 Introduction 

When a homogeneous and isotropic perfect-fluid universe hypothesis is taken, the Ein¬ 
stein equations of general relativity are reduced into a set of linear ordinary differential 
equations, governing the Hubble parameter in terms of a scale factor, known as the Fried¬ 
mann equations, which are of basic importance in evolutionary cosmology. In the recent 
systematic works BEl. we explored a link between the Chebyshev theorem 011] and the 
integrability of the Friedmann equations and obtained a wealth of new explicit solutions 
when the barotropic equation of state is either linear or nonlinear, in both cosmic and 
conformal times, and for flat and non-flat spatial sections with and without a cosmological 
constant. The purpose of the present paper is to present several analytic methods which 
may be used to obtain either exact expressions or insightful knowledge of the solutions 
of the Friedmann equations beyond the reach of the Chebyshev theorem. These meth¬ 
ods may be divided into three categories: roulettes, explicit integrations, and analytic 
approximat ions. 

The idea of roulettes offers opportunities for obtaining some special solutions. Recall 
that many textbooks mention the fact that the graph of the scale factor a{t) against 
time t for a closed Friedmann-Lemaitre cosmology supported by a fluid whose pressure 
is negligible follows is cycloid, thus providing an attractive image of the wheel of time 
steadily rolling from what Fred Hoyle dubbed a Big Bang to a Big Crunch reminding one 
of Huyghens’ isochronous pendulum and inviting speculations, such as those of Richard 
Chace Tolman about cosmologies which are cyclic or cyclic up to a rescaling. More 
recently such models have been called Ekpyrotic following an old tradition of the Stoics. 

Few of those textbooks give details or a derivation of the roulette construction of a 
cycloid and none, as far as we are aware, explain whether or not it is a mathematical fluke 
or whether any plot of the scale factor against time may be represented as a roulette, that 
is as the locus of a point on or inside a curve which rolls without slipping along a straight 
line. 

Evidence that this may be so is provided by what is often called a Tolman Universe, 
that is one in which the fluid is a radiation gas whose energy density p and Pressure P are 
related by p = |P. In this case, for a closed universes, a{t) is a semi-circle which of course 
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may be obtained as a sort of limiting case by “rolling without slipping” the straight line 
interval given by it radius. Further evidence is provided by De-Sitter spacetime which 
contains just Dark Energy. In its closed form, the scale factor has the form of a catenary 
which is the locus of the focus of a parabola rolled without slipping along its directrix. 
In the light of this lacuna in the literature, it is here proposed to provide a general 
discussion, showing among other things that a roulette construction is always possible 
and how in principle to obtain the form of the rolling curve. In order to do so the article 
commences by recalling some of the classical theory of plane curves, mainly developed 
in the and 18**^ century, although its roots go back to the Greeks. As the theory is 
then applied to the Friedmann and Raychaudhuri equations familiar in cosmology it will 
become apparent that the ideas have a much wider relevance. In particular they extend 
immediately to the motion of light, or sound rays governed by Snell’s Law moving in a 
vertical stratihed medium with variable refractive index, and to the familiar central orbit 
problems of mechanics. More surprisingly perhaps we may, as Delaunay hrst pointed 
out in the 19*^^ century, we can also use the same description to describe the shapes of 
axisymmetric soap dims 1316]. In this way we obtain a framework for understanding in a 
unihed geometric fashion a wide variety of iconic problems in physics. 

Recall also that the Chebyshev theorem applies only to integrals of binomial differen¬ 
tials IB El- In cosmology, one frequently encounters models which cannot be converted 
into such integrals. Hence it will be useful to present some methods which may be used 
to integrate the Friedmann equations whose integrations involve non-binomial differen¬ 
tials. Here we illustrate our methods by integrating the model when the equation of state 
is quadratic iniHiE] and the Randall-Sundrum H universe [ini El [13] with a non¬ 
vanishing cosmological constant for which both the energy density and its quadratic power 
are contributing to the right-hand side of the Friedmann equation. A quadratic equation 
of state introduces a larger degree of freedom for the choice of parameters realized as the 
coefficients of the quadratic function. In the flat-space case we show that for the solutions 
of cosmological interest to exist so that the scale factor evolves from zero to inhnity the 
coefficients must satisfy a necessary and sufficient condition that conhnes the ranges of the 
parameters. We also show that when the quadratic term is present in the equation of state 
the scale factor cannot vanish at finite time. In other words, in this situation, it requires 
an inhnite past duration for the scale factor to vanish. As another example that cannot be 
dealt with by the Chebyshev theorem, we study the Randall-Sundrum H universe when 
the space is flat and u-dimensional and the cosmological constant A is arbitrary. There 
are two cases of interest: (i) the equation of state is linear, P = Ap, A > — l,and (ii) the 
equation of state is that of the Chaplygin fluid type, P = Ap— A > —1, B > 0. In the 
case (i) the big-bang solutions for A = 0 and A > 0 are similar as in the classical situation 
so that the scale factor either enjoys a power-function growth law or an exponential growth 
law according to whether A = 0 or A > 0. When A < 0, however, we unveil the fact that 
the solution has only a finite lifespan when A > — 1 and we determine the lifespan 
explicitly. When —1 < A < —1-1-^, the solution is periodic, as in the classical case. In the 
case (ii) with A = 0, we hnd the explicit big-bang solution by integration and display its 
exponential growth law in terms of various physical parameters in the model. Finally, we 
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carry out a study of the big-bang solution for the Friedmann equation when the equation 
of state is of an extended Chaplygin fluid form [TU [151 HH], P = f{p) - Er=i where 
/(p) is an analytic function and Hi,> 0, ai,..., am > 0, for which an integration 
by whatsoever means would be impossible in general. We shall identify an explicit range 
of A for which the scale factor grows exponentially fast and deduce a universal formula 
for the associated exponential growth rate. This formula covers all the explicitly known 
concrete cases. 

The rest of this article is organized as follows. In §2 we review the Friedmann equation. 
In §3 we demonstrate a relation between roulettes and the Friedmann equation and show 
how to use this relation to hnd some special solutions. In §4 we consider the integration 
of the Friedmann equation when the equation of state is quadratic. In §5 we solve the 
Friedmann equation for the Randall-Sundrum II universe. In §6 we establish a universal 
exponential growth law for the extended Chaplygin fluid universe. In §7 we summarize 
our results. In Appendix (§8) we collect some relevant concepts and facts used in this 
article. 


2 Friedmann’s equations 


For generality we shall consider an {n + l)-dimensional homogeneous and isotropic space- 
time with the line element 

= gfj,ydx^dx^ = —dt^ + a^{t)gijdx^dx^, j = 1,..., u, (2.1) 


where gtj is the metric of an u-dimensional Riemannian manifold M of constant scalar 
curvature characterized by an indicator, k = —1, 0,1, so that M is an n-hyperboloid, the 
flat space M”, or an n-sphere, with the respective metric 

gijdx^dx^ = ^ dr^ -F (2.2) 

where r > 0 is the radial variable and df2^_^ denotes the canonical metric of the unit 
sphere in M”, and t is referred to as the cosmological (or cosmic) time. The Einstein 
equations are 

Gfiu + ^g^lu = (2-3) 

where G^y is the Einstein tensor, Gn the universal gravitational constant in n dimensions, 
and A the cosmological constant, the speed of light is set to unity, and is the energy- 
momentum tensor. In the case of an ideal cosmological fluid, is given by 


= diag{-pm. Pm,---, Pm}, (2.4) 

with pm and pm the f-dependent matter energy density and pressure. Thus, assuming an 
ideal-fluid homogeneous and isotropic universe, the Einstein equations are reduced into 
the Friedmann equations 


H 


WnGn k 
n{n — 1)^ a^’ 

- tKP + P) + —, 

n — 1 


(2.5) 

( 2 . 6 ) 
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in which H = ^ denotes the usual Hubble ‘constant’ with / = and p, P are the 
effective energy density and pressure related to pm, Pm through 



(2.7) 


Besides, in view of (12.ip and (12.dh and (12.Th . we see that the energy-conservation law 
= 0 reads 


Pm P 'IT'i^Pm P Pm)H — 0. 


( 2 . 8 ) 


It is readily checked that fl2.5p and fl2.8p imply fl2.6p . Hence it suffices to consider fl2.5p 
and fl2.8p . 


In the rest of the paper, we omit the subscript m in the energy density and pressure 
when there is no risk of confusion. 

3 Friedmann’s equations and roulettes 

In |2], numerous Friedmann type equations arising in a wide variety of applied areas of 
physics, other than cosmology, are integrated as well in view of the Chebyshev theorem, 
which include the equation for light refraction in a horizontally stratified medium, equa¬ 
tions for soap films and glaciated valleys, equation of catenary of equal strength, the 
elastica of Bernoulli and the capillary curves, central orbit equations, equations of capil¬ 
lary curves, equations for spherically symmetric lenses, etc. These equations may also be 
recast in the forms of roulettes as seen below. 

3.1 Rolling without slipping 

In what follows, extensive use has been made of p, Suppose a closed convex curve 7 
enclosing a domain D with boundary 7 = dD lies in a plane H. Physically we think of 
D as a rigid lamina which rolls without slipping along the a straight line L in another 
plane H'. We may suppose the line L to be the a:-axis of a Cartesian coordinate system 
{x,y) for the plane H'. The locus in H' of a fixed point O E D defines a (single-valued) 
bounded (and hence no-mo notonic) periodic graph 7 ' : y = y{x) over L . If 7 is smooth 
and O lies in the interior of D, then 7 will be a smooth simple curve. If O G dD, then 7 ' 
will have cusps where O meets the x-axis. If, by some contrivance, O lies outside D, then 
the curve 7 ' will no longer be a single-valued graph and it will have self-intersections. If 
the point O lies on the curve, one refers 7 ' to a cycloid. Otherwise to a trochoid. Using 
O as origin we may specify the curve 7 by its equation r = r(t). We define r = |r| and 
p = r ■ n , where n is the outward normal of of 7 . Thus r is distance from O to the point 
7 (t) and p the perpendicular distance from O to the tangent of the curve 7 at the point 
7 (t). The relation {p ,r) G is called the pedal equation of 7 with respect to the origin 
O. We also have a relation [y , y') G associated to the graph 7 , where y' = 
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Simple geometry gives 


(p,r) 





If (r, 6) are plane polar coordinates for the plane IT with origin O, then 


P = 





(3.1) 


(3.2) 


Thns given a pedal curve we get a relation {y, y') and if we are fortunate we may be 
able to solve for y{x). Conversely given the graph y{x) we may determine the relation 
(p, r) and if we are fortunate we may be able to solve for p = f{r) and hence determine 
the equation of the pedal curve in in polar coordinates by quadratures. There remains 
the problem of determining whether the pedal curve is among the known plane curves 
catalogued in various books, and if so what are its properties. 

We conclude this section by recalling Steiner’s theorem, which states that 

• the arc length along the rolling curve equals the arc length along the line on which 
it rolls. 


The area under the roulette is twice the area swept out by the rolling curve. 


3.2 The Friedmann equation 

It is convenient in what follows to use units in which SttG^ = 3 and c = 1, and to rewrite 
the Friedmann equation as 

1 + (a)^ = p{a) + (1 — k)a^. (3.3) 

Now we may rescale a and t in fl3.3p by the same constant factor / and transform the 
equation to the form 

1 + {af = fa^pifa) + /^(l - k)a^. (3.4) 

It follows from fl3.ip that the pedal equation of the roulette we seek is given by 

= P^fpifp) + /^(l - fc)F^ (3.5) 

where / may be chosen as we wish. Setting k = 1 and / = 1 gives our basic eqation: 


= p'^pip) ^ 

^ pip) = — 


P‘1 


which may be read in two directions. Either we know p(a), and we deduce the pedal 
equation, or given, a pedal equation, we deduce p(a). In the sequel we shall proceed in 
both directions. 
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Before doing so we recall that in terms of conformal time rj, defined by dt = adr], the 
Friedmann equation may be written (if A: = 1) as 

^ + ka^ = a'^p{a), (3.7) 

which is a hrst integral of a non-linear oscillator equation obtained by differentiating fl3.7p 
with respect to rj, which is equivalent to the familiar Raychaudhuri equation. 

It is readily seen that the integrable cases in view of the Chebyshev theorem include 
the model 

p{a) = aa~‘^ + 13 , (3.8) 

where a, (3, a are arbitrary constants. 

If we introduce [18] the inverse scale factor, or redshift factor by 

6=1, (3.9) 

a 

we hnd 

3 



3.2.1 Steiner’s Theorem 

With this preparation we note that Steiner’s theorem implies that if we are considering 
the scale factor as a function of cosmic time, a = a{t), the arc-length along the rolling 
curve equal to cosmic time and the area J a{t)dt swept out is twice the area swept out 
on the rolling curve. The area under the graph of the scale factor plotted against time 
has no obvious meaning but if instead we consider scale factor as a function of conformal 
time, a = a{p), which we represented as a roulette, then arc-length of the rolling curve 
would be equal to cosmic time and the area under the graph J a{p)dp swept out would 
be equal to cosmic time and as a consequence this would be twice the area swept out by 
the rolling curve. 


3.3 A single-component fluid in a closed universe 

In the simplest barotropic case with k = 1 and 7 7 ^ | we may choose the scaling factor / 
so that 

P-A (3-11) 

In fact John Barrow has pointed out [19] that in this case, if one dehnes 

37-2 

9 = a—, (3.12) 


then the Friedmann equation becomes 

/ 37-2 

dp) V 2 


2 


37-2 


2 


2 


(3.13) 






whence we obtain the Raychaudhuri equation 


/ 37-2 
dr]"^ \ 2 



(3.14) 


and so a general solution for a as a function of conformal time rj is available and cosmic 
time t may then be obtained in principle by quadratures. 

Not surprisingly, in view of the previous subsection and the above discussion, the same 
reduction may be carried out in the context of the slightly more general model fl3.8p with 
a = —87 and 7 7 ^ which leads to the updated equation 


^1 +(fc 

dr,' 


a) 


37-2 




37-2 


(3.15) 


and may be solved as before. 

However we wish to construct the solution as a roulette. The pedal equation is seen 
to be 

r2=/-37^ (3.16) 

which is the pedal equation of the curve 


37-2 
j ' 4-37 


cos 


/ 37-2 

V4-37 



(3.17) 


Note that the 7 = | is a critical case separating solutions which oscillate in conformal 
time and those that blow up exponentially. Another special case is 7 = g, radiation. Both 
require special treatment. 

Some special cases are as follows. 


• For dust, 7 = 1, P = 0, we have 


r = cos^ , (3.18) 

which is the polar equation of a circle of radius where the origin is on its circum¬ 
ference. 

• If 7 = ^, P = we have 

r = cos(20), (3.19) 

which is the polar equation of the Lemniscate of Bernoulli whose Cartesian equation 

is 

{x^ + . (3.20) 

• For radiation, 7 = |, P = d, we have 

r = l, (3.21) 
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which is which is the polar equation of a circle of radius 1 , where the origin is at its 
centre. However the locus of such a rolling curve would be a straight horizontal line. 
This is not consistent with the second-order Raychaudhuri equation (I3.14p whose 
solution is 

a = sin? 7 , t = l — cosr]. (3.22) 

• For stiff matter, 7 = 2 , P = p, we have 

r2cos(20) = l, (3.23) 

which is the rectangular hyperbola 

x^-y^ = l, (3.24) 

with the origin at its centre {x, y) = ( 0 , 0 ). 

• If 7 = |, P = —|p, we have 

r =-(1-|-cos0), (3.25) 

which is the polar equation of a cardioid with origin at the cusp 6 = 71 . 

• If 7 = |, P = —we have 

(3.26) 

which is the polar equation of Cayley’s sextic whose Cartesian equation is 

(4(a;2 + y^)-xf = 27 {x^ + ^ (3.27) 

• If 7 = |, P = —|p, we have 

(3.28) 

which is the polar equation of Tschirhausen’s cubic, whose Cartesian equation is 

27y^ = {1- x){x + 8f. (3.29) 

Tschirhausen’s cubic is also known as L’Hopital’s cubic or Catalan’s trisectrix. 

• For a cosmological term, 7 = 0, 

2 

- = l + cos0, (3.30) 

r 

which is the polar equation of a parabola of semi-latus rectum 2 and origin at its 
focus. 
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We now turn to the special case 7 = |, that is P = — To begin with we assume 
k = 1. Now from fl3.6l) or fl3.3l) po > 1, the pedal equation becomes 


2 2 
r =p — 

sm a 

for some constant a which in polar coordinates becomes 

r = aoe“*“^ 


(3.31) 


(3.32) 


where Oq is a constant. This is a logarithmic spiral unless a = |, in which case it is a 
circle of radius oq. Since Friedmann’s equation reads in this case 


(a)^ = cot^a, (3.33) 

we see that the scale factor increases or decreases linearly with time unless a = | in which 
it is constant. 

3.4 A two-component fluid in a closed universe 

For our next example, consider the locus of a point at a distance A from the centre of a 
circle of radius R. One has 


t = Rt] — Asinr], 
a = R — Acosr]. 


Note that dt = adr] and thus rj is conformal time. Moreover 


da . R sin p 
dt R — A cos 1 ] 


One thus has 

1 + 

The pedal equation works out to be 


42 - p2 + 2Ra 


- R^ + 2Rp, 


(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 


which may also be verified by a simple geometric argument. Comparing fl3.37p with fl3.6l) 
we hnd that 


P = 


m 


A^-R^ 


(3.39) 

a'^ a‘* 

which, if 4 > P is a mixture of dust and radiation with positive energy density. The model 
starts out from a big bang, a increases to a maximum, and the model ends with a big 
crunch. This feature may not be so obvious since the time variable transformation f|3.34p 
is not invertible when A > R. However, in [2], we have carried out a systematic study 
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which shows that the Friedmann equation fl3.37l) is integrable, in view of the Chebyshev 
theorem, to yield a big-bang type periodic solution in cosmic time t. 

Geometrically we have a prolate cycloid. \i A < R, we have a curtate cycloid. The 
energy density of the radiation is negative and this gives rise to a bounce and hence 
oscillatory behaviour. This periodic behaviour is clearly seen in (I3.35p by virtue of the 
invertible transformation fl3.34p when R > A but was untouched in |2]. 

A more general example is given by the pedal equation of an epi- or hypoycloid 


with 


(A+ 25)2 
45(A + 5)’ 


whose parametric equation in Cartesian coordinates is 

x{(3) = (A + 5) cos/3 — 5 cos I 
y{(3) = (A + 5) sin/5 — 5sin ^ 
Comparing with fl3.6p . we see that 


A + B 
B 



A + B 
B 



(3.40) 

(3.41) 


(3.42) 

(3.43) 


aV + (fc - l)a^ = ^ (3.44) 

whence we have 

The last term in fl3.45p corresponds to radiation with a positive pressure, 7 = |- The hrst 
term to a fluid with 7 = §• For particular ratios A/B and hence by fl3.4ip certain values 
of CC we have the following special cases 


• C = |, the Cardoid, 

• C = |, the Nephroid, 

• C = — 1, the Deltoid, 

• C = — 1, the Astroid, 


A/B = 1. 
A/B = 2. 
A/B = -3. 
A/B = -4. 


If fc = 0 or A; = 1 , the first two cases correspond to a fluid with 7 = § with positive 
energy density and negative pressure, the other cases to negative energy density and 
positive pressure. 

Some other cases are as follows. 
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2 ^2 o2 

= —r- h A - 5 

p2 


corresponds to the hyperbola 


x'^ y'^ 


= 1 , 


(3.46) 

(3.47) 


with origin at its centre. \i k = 1 this has stiff matter 7 = 2 with positive energy 
density and radiation whose energy density will be positive if and negative 

if A^ < 52. 


corresponds to the ellipse 


2 A^B^ . 2 

=-^ + A^ + B^ 

pZ 


x^ y‘^ 

+ ’ 


(3.48) 


(3.49) 


with origin at its centre, li k = 1 this has stiff matter 7 = 2 with negative energy 
density, in other words, a phantom, and radiation with positive energy density. 






p ^ 


This is an ellipse with origin at a focus. 

x"^ 

I 2 + 52 

If A; = 1, the energy density is 

P = 


2 9 

^ y ^ 


1 


{2AB)^ (1 + a2)2 

The curve pursued by the scale factor is called an elliptic catenary. 


(3.50) 


(3.51) 


(3.52) 


2 I'AA^B^ ^ ^2 

= I - 5 -4712 

pZ 


-2 


This is an hyperbola but now with origin at a focus. 


A^ 


If /c = 1, the energy density is 


^ (2715)4 (l-a2)2- 

The curve pursued by the scale factor is called a hyperbolic catenary. 


(3.53) 


(3.54) 


(3.55) 
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3.5 ACDM cosmology 


The presently favoured, so-called concordance model has /c = 0 and 

1 A 

corresponding to dust which is known to be integrable. 

The pedal equation is therefore 


r = P\Ip^ 


po 


A 

+ w + 


p^ 


= \Ip+^p^ + p^- 


(3.56) 


(3.57) 

(3.58) 


There is a simple formula for the scale factor as a function of time but hnding the 
curve of which it is the roulette appears to be difficult since one has to solve a quartic 
curve for p. 

The reader will have noticed that many of the explicit examples of roulettes giving 
rise to solutions of Friedmann’s equations descibed above have involved “exotic equations 
of state” with components having for example negative energy densities and/or pressures 
exceeding energy densities in magnitude. These were often excluded from traditional 
texts on cosmology on the grounds of lack of physical interest. However with the growing 
observational support for the presence of dark energy the literature has become a great 
deal less inhibited (see e.g. ra), and we have therefore felt justihed in including as many 
explicit cases of the roulette construction as we could hud. 


4 Quadratic equation of state 

It is well known that the flat-universe {k = 0) Friedmann equation in the linear equation 
of state case 

P = Ap, (4.1) 

where A is a constant satisfying the non-phantom condition A > — 1, may be integrated for 
arbitrary values of the cosmological constant A to allow big-bang solutions with a(0) = 0, 
so that when A = 0 the scale factor a{t) grows following a power law, when A > 0, an 
exponential law (a dark energy regime), and, when A < 0 the scale factor a{t) oscillates 
periodically. 

In this section, we consider the more general situation [3 Eli when the equation of 
state follows the quadratic law 

P = Po + Ap + Bp^, (4.2) 

which may be viewed as a second-order truncation approximation of the general equation 
of state, P = P{p). We note that cosmologies governed by quadratic and some other 
more extended forms of equations of state were investigated earlier in [211 [22l [23]. For 
convenience, we again assume non-phantom condition, A > —1, throughout our study. 
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4.1 Restrictions to parameters in non-degenerate situations 

Inserting (14.2^ into the law of energy conservation, 


p + n(p + P)- = 0, 
a 


we obtain 


J = 


dp 


= ln(Ca-”), 


(4.3) 

(4.4) 


Fo + (1 + ^)p + Bp‘^ 

where C > 0 is a constant. 

For convenience, use D to denote the discriminant of the quadratic 

Q{p) = -Pq + (1 + ^)p + Bp^. (4.5) 

That is, D = {l + Af - AEP^. 

For F 7 ^ 0, there are three subcases to consider. 

(i) F = 0. Then J = — where po = — Hence, absorbing a possible 

integrating constant, we have 

1 


B{p - Po) 


= ln(C'a-’^). 


(4.6) 


We are interested in a scale factor a initiating from sufficiently small values and 
evolving into arbitrarily large values (‘big bang cosmology’). If F > 0, then —p > 0. 
Using this fact and p > 0 in fl4.6l) . we see that small values of a are not allowed. If 
F < 0, then fl4.6p assumes the form 




= ln(Ca-”). 


(4.7) 


To allow a to be small, we obtain the restriction |F|p—4(1 + 4) > 0, which prohibits 
a to assume large values. In conclusion, the case F = 0 is ruled out. 


(ii) F > 0. Then Q{p) has two roots: 


Pi,2 — 


-(l + 4)±^/F 


2F 


(4.8) 


resulting in J = F 2 In 

° P-P2 


Hence 

Pi ± Ca~'^^p2 

P = 


(4.9) 


where F > 0 is an updated constant. Inserting the epoch a = 0 into (14.bh . we get 


, 1 + 4 + \/F 

p|a =0 = Po = P 2 = -^Tgl- >0, F < 0. 


(4.10) 
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In other words we find a sign restriction for B. Inserting a(cx3) = cxo into fl4.9p . we 
get 

(4.11) 

(4.12) 


p(oo) — Pi — + ^] ~ ^[l + + A\B\Po) > 0. 


Therefore we must have 

Fo<0. 

Using Pi < p{t) < p 2 (Vf), we see that fl4.9p is made precise: 

_ pia"'^ + Cp2 
P g^nVD (J 

(iii) F < 0. Then Q{p) has no real roots. We can integrate fl4.4p to get 


2B \B 

arctan 






2p + BA I = i„(ca-"). 


(4.13) 


(4.14) 


Since the right-hand side of fl4.14p stays bounded, we see that small values of a are 
not allowed. Thus the case F < 0 is also ruled out. 


4.2 Global solutions 

Hence we now aim at integrating the fiat-space {k = 0) Friedmann equation 




:P + 


2A 


n{n — 1) n{n — 1) 
in the case D > 0 only. For this purpose, we may insert fl4.13p into fl4.15p to get 

2 r _ ny/D 


(4.15) 


L^gnVlJ + L2 


where 


Li = 


n{n — 1) 

Thus, an integration of 04.161) gives us 


2 2C 

(8vrG„pi -l- A), L 2 = - -T{9>TiGnp2 + tV). 


n(n — 1) 


t = 1 = 


da 


Lio^^ + L2 a 


“ + ^ (u = a"A, 


ny/D J V Liu + L 2 u 

for which the Chebyshev theorem is not applicable in general. 


(4.16) 


(4.17) 


(4.18) 
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If Li = 0, then L 2 7 ^ 0 otherwise pi = p 2 and D = 0 which is false. Thus L 2 > Li = 0. 
In this case we can set u = L 2 V^ — C to get 


t 


2 f L2V^ 
nVDj L2V^-C ^ 



+ Cl, 


(4.19) 


where Ci is an integration constant and 




V = 


(4.20) 


In other words, fl4.19p and fl4.20p give us the general solution of the Friedmann equation 
when Li = 0 . 

From this solution we see that the initial condition a = 0, or 


C 


V = 





(4.21) 


cannot be achieved at any hnite time but can only be realized ai t = — 00 . More precisely, 
from (14.191) and (I4.2np . we may deduce the property 


a{t) = O eV 1 ? 


t —)■ —cxo. 


(4.22) 


Furthermore, from fl4.19p and fl4.20p . we may also deduce the asymptotic behavior 

a{t) = O t —)■ cxD. (4.23) 


Now assume Li 7 ^ 0. From fl4.16l) we have Li, L 2 > 0. Hence Li > 0. Thus L 2 > Li > 
0. In this case we can set 


V = 


Liu + L2 

Consequently, we obtain 

2{CLi-L2) 


u + C L 2 V^ - C 

or u = 




- 1 


t = 1 = 

1 


n 


y/D 


{Liv^ - 1){L2V^ - C) 


L 2 

Liv^ - 1 Ti 


dn 


(4.24) 


n 


y/D \ y/TCi 


In 


V + 


y/Li 


VLt 


+vl"' 




»’ + >/£ 


+ Cl, 


(4.25) 
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where Ci is an integration constant. Since 


V = 


+ L 2 


(4.26) 


we see that a = 0 corresponds to v = J ^ like before. 


From fl4.17p . we have 


C n{n — 1) 


(87rG„p2 + A) > 


n{n — 1) 


(SttG^pi + A) — Ai, 


(4.27) 


which leads to f Consequently we see in view of fl4.25p and fl4.26p that the scale 

factor a cannot vanish at any finite time but at t = — 00 . More precisely, we have 


a{t) = O ( eV c 


I 


t —>■ — 00 , 


(4.28) 


which is the same as fl4.22p and may hardly be surprising. Similarly, we have 

a{t) = O , t ^ 00 . 


(4.29) 


This is an important conclusion since we have achieved an exponential asymptotic growth 
law. It is worth noting that the growth rate y/Li is given by the explicit formula 


Li = 


47rGr, 


n{n — 1) \ \B\ 
B <0, Po < 0. 


(1 + A)-V(1 + A)2-4|PPo|1 + A), 


(4.30) 


In summary, we can conclude: 

(i) A nontrivial quadratic equation of state given in fl4.2p dictates the necessary and 
sufficient condition 

Pq ^ 0) P < 0, (1 + A)^ + 4|P|Pq > 0, 

(4.31) 

^ ((1 + A)- y(l + Af - 4|SF„|) + A > 0. 

in order to allow the scale factor a to evolve from the epoch a = 0 to the epoch 
a = 00 . 

(ii) The quadratic equation of state fl4.2|) may never allow a hnite initial time at which 
the scale factor vanishes. 
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(iii) In the critical situation 


((1 + A) - ^{l + Ay-A\BPo\) + A = 0, (4.32) 

the scale factor a grows following a power law. On the other hand, in the non-critical 
situation 

((1 + A) - V(1 + A)2-4|5Po|) + a > 0, (4.33) 

however, the scale factor a obeys the exponential growth law (I4.29p ~ (l4.3np . In par¬ 
ticular, in this latter situation, the cosmological constant A is permitted to assume 
a suitable negative value when Pq < 0 . 


We may compare our results with the a-fluid model studied in PI for which the 
equation of state is 

P = —(a -|- 1 )pa P (^P — (a -|- 1) —, (4.34) 

Pp 

where a > 0 is a parameter, pp is the Planck density, and p\ the cosmological density. It 
is seen that fl4.34p is consistent with a part of the necessary and sufficient condition just 
derived. Furthermore, the positivity condition for the discriminant of fl4.5p . i.e., P > 0, 
translates itself into the form 


/p i 

Pp 4’ 


(4.35) 


which does not involve a, in particular. Since pA is of order 10 and pp of order 10®® 
(cf. [2l]), the condition (I4.35p is well observed too. 


4.3 Degenerate situation 

In the degenerate situation, P = 0, (14. 4 p gives us the relation 

|Po + (1 + A)p\ = (4.36) 

Since we are interested in the case when a may assume arbitrarily large values as / —)■ cxo 
and p > 0, we again arrive at the necessary condition 

Po < 0. (4.37) 

Hence 

P = + |Po|) . (4.38) 

In other words, the presence of Pq simply and unsurprisingly adds a shift to the cosmo¬ 
logical constant. 
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4.4 Polytropic equation of state case 


It will also be interesting to extend our study in this section into a slightly more general 
situation where the equation of state is of a polytropic one given by |25l ES] 

P = Ap + Kp'^, A >—1, 7 7^1, 7 > 0, (4.39) 

where k ^ 0 otherwise the equation of state becomes the well-studied linear case. 
Inserting fl4.39p into the equation of energy conservation, we have 

(1 + A)p^-^ + K = , (4.40) 

where C > 0 is an integration constant. 


7 < 1. Then fl4.40l) indicates that the epoch a = 0 gives rise to the initial value 
Po = oo and the epoch a = oo gives rise to the condition k < 0 and the limiting 
value 


Poo 



(4.41) 


(ii) 7 > 1. Then fl4.40p implies that the epoch a = 0 gives rise to the condition k < 0 
again and the initial value 

and the epoch a = oo gives rise to the limiting value poo = 0. 


Note that the necessary condition k < 0 is analogous to the necessary condition B < 0 
stated in fl4.3ip for the quadratic equation of state problem. 

Applying (I4.4np in the flat-space Friedmann equation with A = 0 for simplicity, we 
arrive at the integral 

1. 

4 ( ' (1 + A)-^ t = I= [ a-^ + \k\)~^ da. (4.43) 

\n{n-l)J J 

In view of the Chebyshev theorem, we know that the integral I is elementary when A and 
7 are arbitrary rational numbers. Here we omit presenting the solutions. 


5 The Randall—Sundrum II cosmology 


In this section we consider the braneworld cosmological Friedmann equation ITO 

k 

H—r = otip + (y^2P^ A 


(5.1) 


in the Randall-Sundrum II universe mi [12113] where 


IbTrGn 


CKl = 


/ T \ 1 

n[n — 1) y 


(IQnGn+i 


n 


X = 


2 A 


n{n — 1) 


(5.2) 


are specihc physical parameters. 
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5.1 Linear equation of state case 

In view of the linear equation of state 


P = Ap 


(5.3) 


and the law of energy conservation fl4.3p . we have 

p = Ca-^^^+^\ 1 + A>0, C>0. (5.4) 

Inserting fl5.4p into fl5.ip and restricting to the flat-space case k = 0, we easily obtain 
the integration 


t Cl — / — 


da 


J /3ia-A^+^) + y 

^^a2n(l+A) + p^aA^+A) + ^2 


dti 


n{l + A)J ^Xu'^ + + {32 


(u = 


= 


n{l+A)pi 

1 

n(l+A)%/X 


A = 0, 


^ In + ^a^Ai+A) + §Aan{i+A) + ^ y > g, 


n(l+A)y/\P\ 


arctan 


tCi+ai Pi 
2|A| 


y_a2n(l + A) + ^(_3jQ,n(l + A)+^2) J ’ 


A < 0, 


(5.5) 


since it is of the Chebyshev form, where Pi = Cai and P 2 = C‘^a 2 and Ci is an integration 
constant which may always be chosen to allow the big-bang initial condition a(0) = 0. 
Indeed, from fl5.5p . we have 


6*1 = 


n(l+A)P2 ' 

1 .\ii I §1 ^ Si] 

2A ^ A ’ 


n(l+A)\C 

_ 1 

n(l-|-A)y^ 


2A ' V A 

arctan 


A 


2v1 ^J ’ 


A = 0, 
A > 0, 

A < 0. 


(5.6) 


From the above results we deduce, as t —>■ 00 , the growth laws 


a(t) = O ) , A = 0; a(t) = O (e 


,y/Xt 


A > 0. 


(5.7) 


These two situations are analogous to those in the classical situation mm- 

However, it may be a surprise to see, when A < 0, the solution will be of a finite 
lifespan or periodic, instead, depending on how far the parameter A is away from the 
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phantom divide line [28l |29l [30l |3T1 [32], A = —1. To see this fact more transparently, we 
rewrite the Friedmann eqnation as 


(^) =n‘^{l +Af\X\{ui-u){u + U 2 ), 


where 


Ml = 


2|A| 


and u = Using m( 0) = 0, we solve 

du 


to get 


dt 


at = arctan 


= a^y (mi - m)(m + M 2 ), a = n{l + A)^/\^, t > 0, 


I M-|(Mi-M2) 

W(Mi -m)(m + M2) 


+ arctan 


1 (mi - M 2 ) 


t > 0. 


2 y/UlU2 

The above relation remains valid nntil t = to > 0 (say) when u = Ui. Thns 


to — 


a V 2 


TT 


+ arctan 


1 {ui -U 2 ) 

2 


(5.8) 


/?i + \lI3l + 4/?2|A|^ , M2 — -^1^1 f3l + 4/?2|A| — /3i^ , (5.9) 


(5.10) 


(5.11) 


(5.12) 


Beyond to the eqnation (15.101) is invalid. Instead, m reaches a ‘stagnation’ point, where 
M(to) = 0 and u(to) = Ui, and starts to descend. Hence it must evolve through another 
branch of the Friedmann equation. 


— = -a\/ (mi - m)(m + M 2 ), a = n{l + A)^/\^, t > to, (5.13) 

dt 

which may be integrated similarly. Actually, since both fl5.10p and fl5.13l) are autonomous, 
the solution of fl5.13p satisfying u(to) = ui may simply be obtained from that of (15.101) 
through replacing the variable t with t G [0,to] for the solution of (15.101) by 2to — t with 
t G [to,2to]. In particular, u{2to) = 0. 


5.2 Solutions of finite lifespans 

We now show that the solution ceases to exist beyond 2to when 

A >-! + -. (5.14) 

n 

In fact, from (I5.13p . we have 


lim u(t) = —aJuiU 2 < 0. 


(5.15) 
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From the relation between the scale factor a and the solution u we hnd 


= n ^ (5.16) 

n(l + A) 

so that 

lim ait) = —oo, A > —1+—; lim ait) = —a\/uiU2 <0, A = —1+ — . (5.17) 
t^{2to)- n t^{2to)- n 

On the other hand, if the solution would exist beyond 2to, we should have a{t) > 0 when t 
passes 2to where a = 0. Thus (I5.17|l indicates that ait) would suffer from a discontinuity 
at t = 2to- In other words, the solution to the Friedmann equation cannot exist beyond 
2 to which may well be called the lifespan of the solution. 

If the coupling constant A satishes 

71 <-! + -, (5.18) 

n 

then (15.161) implies that a(t) —)■ 0 as t —?• (2to)~- Thus we may obtain the solution of the 
Friedmann equation over [2to, 4to] by shifting the solution over [0, 2to] obtained above and 
maintain the continuity of ait) at f = 2to- Such a procedure can be repeated so that a 
periodic solution of period 2to is constructed. 

Thus, we conclude that, when the cosmological constant is negative, the solution a of 
the Friedmann equation evolving from the initial condition a(0) =0 has either a hnite 
lifespan T = 2to or is of period T = 2to, according to whether the coupling constant A 
fulhlls the condition fl5.14p or fl5.18p . 

For later convenience, we document the lifespan or period T here explicitly in terms 
of the original physical parameters: 


T 


1 In — 1 
1 + A\I 2n\A\ 


^TT + 2 arctan 


Gn I n ~ 
2Gn+i\l 2(n-l)|A| 


A < 0. 


(5.19) 


Thus we have unveiled a peculiar situation, A > —1 + 4^ when the solution has only a 
hnite lifespan. This situation is in sharp contrast with the classical case with A < 0 that 
the big-bang solution is periodic for any A > — 1 . 

To end this subsection, it will be instructive to check that, subject to the equation of 
state fl5.3p . the phantom world is still spelled out by the condition 1 4 - A < 0 as in the 
classical situation. In fact, for 1-l-A 7 ^ 0, the energy conservation law gives us the relation 
p = as before where C > 0 is a constant. Inserting this into the Friedmann 

equation fl5.ip with /c = 0, A = 0, we have the integration 

t = I= [ . (5.20) 

J aVC'aia-"(i+^) + G^a2a-^^(^+^i 

Hence the scale factor a = a{t) is implicitly given by 

^^n(i+A) a’^(i+^)(0) + a2 + ^n(l -4 A)ai1^ , t > 0. (5.21) 
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Consequently we see that the initial condition a(0) = 0 is permitted when 1 + ^ > 0 so 
that a = a{t) grows following the law 


a{t) = O (t"(i+^) 


t —)■ oo. 


(5.22) 


If 1 + ^4 < 0, however, the initial condition a(0) =0 is not permitted and the Big Rip 
[331 [3^ happens at 


to — 


2 A 


1 A\ai V 

Ic/ 


a"-(i+^)(0) + a2. 


(5.23) 


In particular, A = —1 still corresponds to the phantom divide line |2Sl 123 EDI EB E2] 
indeed as in the classical situation, as anticipated. 

In [35], the following slightly more general Friedmann equation 


— ttp + j3p^ H—j 


(5.24) 


is considered {k = 0 ,A = 0), where a,l3,'y > 0 are constants. Although this equation 
cannot be integrated, some asymptotic analysis establishes that any solution with a(cx 3 ) = 

oo grows according to the law a(t) = O j a.s t ^ oo, as before, subject to the 

equation of state P = Ap (A > —1). 


5.3 Chaplygin fluid case 

Consider the equation of state 

B 

P = Ap -, A>-1, R>0, (5.25) 

P 

of a generalized Chaplyin fluid model [36l [ST] [38]. Inserting fl5.25p into the law of energy 
conservation, fl4.3p . we obtain 

P = " , (5,26) 

where C > 0 is an integration constant. Substituting fl5.26p into the braneworld Fried¬ 
mann equation fl5.1ll with k = 0 and A = 0 and integrating, we arrive at 

t = I^ j (a, + bY+a2 b = (5.27) 

To proceed, we note that (I5.26p gives us 

a = (p^ — 6 )“ 2 'i(i+^) , p^ > b. (5.28) 

Combining (I5.27p and (I5.28p . we have 
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which is not of the Chebyshev type. Nevertheless, set 


u = 


P 


Q!l + Ot2p 


Then 0 < m < and 


P = 


aiu 


1 — a2U^ ’ 

which renders the integral fl5.29p into a rational one, 


I = 


u 


dM, 


+ A) J — bai ^(1 — a2U^y 

which can be integrated separately in the parameter regimes 

(y.\ cx\ (y.\ 

respectively. 

We begin by displaying the result in the first regime as follows: 


I = 


n(l + A)\/h 


02 + 


0^1 

7^ 


ai 


arctanh | ( 0:2 + ^ I u 


- 02 - 


7^ 


arctanh ( 0:2 — 


Oil 

7^ 


(5.30) 


(5.31) 


(5.32) 


(5.33) 


Oil 


u\)+Ci, a2> (5.34) 


where Ci is an integration constant. It is clear that we may choose Ci suitably to allow 
the big-bang initial condition a(0) = 0 or m( 0) = On the other hand, the growth 

condition a(cxD) = 00 gives rise to p{oo) = \/b. Consequently, we arrive at 

«2 + , (5.35) 

as anticipated. In other words, the first term on the right-hand side of fl5.34p dominates 
as t ^ 00 , which leads to the following explicit asymptotic estimate: 



i{t) = O (e”"*) , t ^ 00, 


where 


a 


= Vb { a 2 + 


ai 

Th 


B 


l + A 


IOttG, 


n 2 


n+l 


n 


+ 


WnGn l + A 


n{n — 1) V B 


under the condition 


«! < a2\/b or Gn < IbTT 


B 


l + A 


(5.36) 


(5.37) 


(5.38) 
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Likewise, in the third regime, we have 


Q!in(l + A)I = ^ + 0^2 


0^1 / 0L\ 

“7^ 


arctanh 


arctan 


«! 

\/h 

0^1 


+ 02 U 


y/h 


- Q!2 M + Cl, 


where a 2 < 


0^1 

7^ 


n 


or 167r( 1-- 1G7 i\/ 


5 


< G„. (5.39) 


With fl5.39p . we have the same asymptotic state fl5.35|) and we see that fl5.36p and 
fl5.37p still hold. 

We omit here the critical case \/ha 2 = oii. 

It is interesting to note the different ways the Newton constants and and the 
constants A and B make their contributions to the dark energy, a, as expressed in (I5.37p . 
It will also be enlightening to compare fl5.36p - fl5.37p with the formulas [2] 


cr = 4 


irGr 


B 


n{n — 1)/ \1 + 

obtained for the big-bang solution of the classical Friedmann equation 


= 


IQnGr. 


( -t \Pi 

n[n — 1) 

which may simply be obtained from fl5.37p by setting G„+i = 0. 


(5.40) 


(5.41) 


6 Universal growth law for the extended Chaplygin 
universe 

Recall that the generalized Chaplygin fluid is governed by the equation of state 
P = Ap -, A > —1, R > 0, 0 < a < 1, 


( 6 . 1 ) 


whose cosmological implication, especially its relevance to dark energy problem, has been 
studied extensively. It is direct to see from fl6.ip that there holds the scale factor and 
energy density relation 

( 6 . 2 ) 

^ 1+A ^ ' 

In view of the Chebyshev theorem, we know that the Friedmann equation with k = 0 
and A = 0 is integrable [2] when a is rational. Nevertheless, for any real a, we derived in 
[2] the following explicit universal growth law: 


a{t) ~ Coe”" , t —)■ cx); a = 4 


ttG. 


n{n — 1)/ yl -l- A 


B \ 2 ( 0 + 1 ) 


(6.3) 


26 


























Inserting fl6.3p into fl6.2p . we have 


/ B \ 

, t^oo. (6.4) 

Unsuprisingly, poo is the only positive root of the function 

h{p) = P + p = p +Ap- (6.5) 

Suggested by the above results, we establish in this section a universal growth law in 
much extended situations, without resorting to an integration of the Friedmann equation. 


6.1 Extended Chaplygin fluids 

In [m US] and [IB] , the extended Chaplygin fluid model given by the equation of state 


P 




i=i 


B 


( 6 . 6 ) 


is studied in the specihc situations when Ai = A > 0 {I = 1,...,M) and Ai = j 
{I = 1,...,M), respectively, by numerical means. Motivated by their work, here we 
consider the extended Chaplygin model governed by the following general equation of 
state: 


N 

p = f(p) - E 

m=l 


Bm 

pOLm ^ 


p > 0, 


( 6 . 7 ) 


where i?i,..., Bj^ > 0 and cxi,..., cktv > 0 are constants and / is an analytic function 
satisfying 


/(0)=0; /'(p) >-1, p>0; p + /(p)-^cx) as p cx). (6.8) 

A direct consequence of the above assumption on / is that the function 

^ R 

ft(p) ^ P + /(P) - E ^ (6-9) 

m=l P 

has exactly one positive root, say poo > 0, since h(p) strictly increases in p > 0. 

It is clear that our assumption on /, motivated from the results in [2] and discussed 
above, covers all the examples considered in [IH UH 116] . 

6.2 Universal growth law 

Rewrite the function h in (16.9p in the form 

h{p) = ip-Poo)hi{p)- ( 6 . 10 ) 
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Since h{p) < 0 for p > 0 small, we see that hi{p) > 0 for all p > 0. Inserting this into the 
law of energy conservation fl4.3p . we have 



dr 

hi(r)(r - Poo) 



( 6 . 11 ) 


where po = p(to) > 0 and to is an initial time when a(to) > 0 and p = p(t), a = a{t) with 
t > to- Note that we are dealing with an expanding universe such that po > p > Poo- 
In order to extract information from fl6.1ip . we rewrite the integrand of it as 


hi(r)(r - Poo) 


hi{poo){r - Poo) r-poo\hi{r) hi(poo) 
1 

+ q{r), r > 0, 


hiipoo)ir - poo) 

where q{r) is a smooth function in r > 0. Inserting (16.1211 into (16.lip , we arrive at 

1 


( 6 . 12 ) 


hi(po 


In Ip - Pool = Q{p) + hi (a ") , 


where 


(6.13) 


Q{P) = 


hi (poo) 


ln|po-Pool+nln(a(to)) - / q{r) dr, Poo < P < Po, (6.14) 


'PO 


which is a bounded function. Combining (16.131) and (16.141) . we have 

P = Poo + e^hP oo )Q{p) g^rihi{p 

OO ) 

Inserting (I6.15p into the Friedmann equation (14.151) . we get 


where we set 


= Ao + ^ _ 167^^ 

a J n[n — 1) 


^ 167rG„ , 2A ^ ^ 

Ao — —- ttPoo H-^-7T > 0. 


n{n — 1) n{n — 1) 
In view of (I6.16P and (I6.17p . we obtain the integration 


t = 


da 


ai/ Ao + « e^hPoo)Q(p)Q-nhi(poo) 


da 


\/A()( 


tV"Ao + a\/Ao ( 


da 




In a + R{a), 


(6.15) 


(6.16) 


(6.17) 


(6.18) 
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where 


R{a) = 


aVAo + a e^hPoo)Q(p)(2 nhi{p^) \/Ao a 

1 f ae^Ap oo CXD ) 

2 ' '^y/Ao + a^(a)e'"i(poo)g(p)a-n/ii(p<x.) 


da 


da, e(a)e[0,l], (6.19) 


which is bounded for a near oo since hi(poo) > 0. Consequently, in view of (16.181) and 
fl6.19p . we deduce the universal growth law for the scale factor a as follows: 


a{t) = O ^e^*j , f —)■ oo. 


( 6 . 20 ) 


where Aq is given in fl6.17p . 


6.3 Examples 

It will be enlightening to display a few special examples which allow explicit calculations 
of Poo, and hence, Aq. 


(i) For the classical generalized Chaplygin model fl 6 .ip . we see that poo is given by fl6.4p . 
Hence 

16xG„ f B ^ 2A 

° ?7,(?7, — 1) \1 + H/ ?7,(?7, — 1)’ 

which generalizes the corresponding formula obtained in [ 2 ]. 

(ii) Consider the extended Chaplygin fluid governed by the equation of state 

B 


P — Aip + H 3 P-. 

P 


It is direct to see that 

Poo 

Consequently we have 
167rG„ 


2A 


3 L 


\/(l + Ai)^ + 4 ^ 3 !? — (1 + Ai) 


An — 


2B 


n[ 


+ 


2A 


i{n - 1) (1 + Ai) + v/(l +Ai)2 + 4A3H ) n{n - 1)' 

In the limit A 3 —)■ 0, fl6.24p recovers fl6.2ip for a = 1. 


( 6 . 22 ) 


(6.23) 


(6.24) 
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(iii) As another simple example, let the equation of state be 


P = Ap-Bo-—. 

P 


Then we have 


Poo — 




2(1 +A) 

which allows us to return to fl6.2ip with a = 1 again in the limit i?o —)■ 0. 


(6.25) 


(6.26) 


There are other explicitly solvable cases of interest as well. For example, we consider 
the equation of state 

P = Ap-^-^, a2>ai>0. (6.27) 






To solve the associated equation 


(1 +A)p - — - — = 0, 


(6.28) 


we require 1 + 0^2 = 2(^2 — ai) or 0^2 = 1 + 2ai. Hence ai > 0 is allowed to be arbitrary, 
which contains the case (iii) above as a special example where ai = 0. In other words, 
we hnd an explicit solvable general situation with the equation of state 


P = An- —- 

P pa pl+2a ’ 

With fl6.29p . we obtain p^o (hence Aq) explicitly: 


a > 0. 


(6.29) 


poo 


Bi + ■sj Bf + 4(1 + A)B2 

2(1 +A) 


l + a 


(6.30) 


Similarly, with 


we have 


B 


P — Alp + A 2 P “--, Q; > 0, 

P“ 


poo 


^7(1 + ^ 1 )^ + 4y42-B — (1 + Ai) 

2X 


1 

1 + Q 


(6.31) 


(6.32) 


which covers (ii) with a = 1. 


7 Conclusions 

Our results in this work are summarized as follows. 
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(i) Friedmann equation and roulettes. We have shown that every solution of Friedmann 
equation admits a representation as a roulette, that is the locus of a point in the 
plane of a curve which rolls without slipping on a straight line. A method is given 
for hnding the curve. A well-known example is the scale factor of closed universe 
containing pressure free matter. The scale factor is a cycloid, i.e. the locus of a 
point on the circumference of a circle, which rolls without slipping on a straight line. 
Adding radiation, one hnds that curve is still a circle but the point moves outside 
its circumference. Many other explicit examples are given, some involving exotic 
equations of state which have been of current interest to cosmologists in connection 
with dark energy and the pre-inflationary universe. 

The Friedmann equation, or its equivalent, also arises in other contexts in physics 
not directly connected with cosmology such as central orbit problems and geometric 
optics. Examples of the use of roulettes are in this wider setting are given. 

(ii) Quadratic equation of state. Although the flat-space Friedmann equation is not 
in a binomial form so that the Chebyshev theorem is applicable, a necessary and 
sufficient condition is unveiled for the coefficients of the equation of state to allow 
the scale factor a to evolve from zero to inhnity. A distinguishing characteristic 
is that the presence of the quadratic term in the equation of state prohibits a to 
vanish at a hnite initial time. In other words, the epoch a = 0 may only happen 
at t = —cx). In the general non-critical situations, a grows exponentially even when 
the cosmological constant A assumes a negative value above a specihc critical level. 

(hi) Randall-Sundrum II universe. In the Randall-Sundrum II universe, there are two 
interesting flat-space cases: When the equation of state is linear, it is shown that 
the scale factor a grows following a power law or an exponential law according to 
whether the cosmological constant A is zero or positive, as in the classical situation, 
but when A is negative, however, a new phenomenon occurs so that in an explicitly 
given regime away from the phantom divide line, the big-bang solution satisfying 
a(0) = 0 has only a hnite lifespan, which never happens in the classical situation, 
although near the phantom divide line the solution is periodic, as in the classical 
situation; when the equation of state is that of a Chaplygin huid, we explicitly 
derive the big-bang solution and describe its exponential growth pattern in terms 
of various coupling parameters, when A = 0. We see as in the classical Chaplygin 
huid model case that a small amount of exotic nonlinear matter would lead to a 
large amount of presence of dark energy near the phantom divide line, as realized 
by the exponential growth rate formula. 

(iv) Universal growth formula. For the extended Chaplygin huid universe for which the 
equation of state is decomposed into the sum of a positive analytic function and a 
negative inverse power function, an integration is impossible in general. To tackle 
the problem, an analytic method is introduce which enables us to extract all the 
key parameters that give rise to a universal formula for the exponential growth rate 
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of the scale factor. Such a formula covers all known formulas derived in concrete 
situations. 
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8 Appendix: Pedal equations and their inversion 

If a curve in the plane is given in polar coordinates (r, (p) by 

u = ^ = u{(j)), (8.1) 

and if p is the perpendicular distance from the origin to the tangent at the point (r, (/>), one has 

2 


1 f du\' 2 

^ = wJ • 

Conversely, if we are given the pedal equation of a curve in the form 

4 = 

p^ 


( 8 . 2 ) 


(8.3) 


the problem of finding the polar equation of the curve amounts to solving the differential equation 

2 


1 /'du\‘ 2 


Equation (|8.4I1 arises in many other contexts. Here are two examples. 


(8.4) 


If a particle moves in a spherically symmetric potential H(i) per unit mass and has 
conserved energy per unit mass £ and conserved angular momentum per unit mass h then 
its orbit satisfies (18.4p with 




(8.5) 


The scale factor a(t) of Friedmann-Lemaitre universe with energy density p{a) (including 
any contribution due to a cosmological term) satisfies 


a2 1 


. „ , . (1 — fc) 

2 H— K — 4:TTGp{a) H- K —, 

Z rtZ /^Z 




which may be brought to the form (18.41) with u = ^, cp = dp = adt, and 

P{u) = 4:TrGa‘^p{a) + (1 — k)a‘^. 


( 8 . 6 ) 


(8.7) 
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8.1 Affine and rescaling of angle 

Suppose ui = f{(p) satisfies 

(^) = ( 8 . 8 ) 

Let U 2 = yl/(A(/>) + B. Then 

(^) +^2 = ^ 2 (« 2 ), ( 8 . 9 ) 

with 

P2{u 2) = A^X^Pi (^ - -S) + 2\^ABu2 + (1 - X^)ul - X^A^B^. (8.10) 

We may compose two such tranformations to get a third: 

2I3 = 2I22I1, A3 = A2A1, B 3 = B 2 + AsBi . (8.11) 


Thus the set of central potentials is acted upon by a three-parameter group. The first two for¬ 
mulas are commutative but the second is not. In the above the second transformation A 2 ,B 2 , X 2 
follows the first Ai, i?i, Ai. 

A matrix representation may be given 



/A 0 B\ /u\ 

Uo 


( 8 . 12 ) 


which reveals the three-dimensional group as corresponding to A times the semi-direct prod¬ 
uct X corresponding to A and B. Acting on P{u), we have 

P{n) ^ A^A^ (^P Q + s) - Q + + u\ (8.13) 


Or if Q{u) = P{u) — V?, then 

g(u) ^ A^A^g Q + s) , (8.14) 

which says that Q{u) transforms under pull-back composed with mutiplication by A^A^. 


8.2 Non-linear transformations 

The simplest of these is to raise n to a power: If 



U2 = {uiT^ 


(8.15) 

then 




P2{U) 

= Pi 

-b (1 — n^)u^. 

(8.16) 

The case n = 1 is particularly 

simple: 




P2{U) = U^Pl (" 

;)■ 

(8.17) 
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8.3 Chebyshev’s theorem 

From (18.41) we deduce that 


du 


J ^Piu) - V? ■ 

Chebyshev’s theorem states that for rational numbers p, q,r (r 7 ^ 0) 0 and nonzero real numbers 
a, [3, the integral 

1 = j.- (a + dx = / (ax5 + dx (8.19) 

is elementary if and only if at least one of the quantities 


p+i. P + 1 , 

— q, — —Hg, 

is an integer. In our case q = and so integrability in finite terms is possible if 

P{u) = + au'^ + , 

where is an integer or half integer. 

8.4 Examples 

8.4.1 Circles 

With respect to a point distance A from a circle of radius R 


r'^ = A^ -R^ + 2Rp, \ = 


AR^u 


2nA 


{I + [R^^ - A^)v?f' 


( 8 . 20 ) 


( 8 . 21 ) 


( 8 . 22 ) 


8.4.2 Conics 

An ellipse 

with respect to the focus 

and with respect to the centre 

A hyperbola 

with respect to the focus 

in this section should not be confused for the perpendicular distance from the origin to the tangent. 


A2 + 52 A 

(8.23) 

p^ B'^ (y2Au ’ 

(8.24) 

1 A‘^ + B^ 1 

(8.25) 

A2 52 “ ’ 

(8.26) 

' = Vi M 

p 2 52 2Au)' 

(8.27) 
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and with respect to the centre 


(8.28) 


A parabola 

with respect to the focus 


with rsepect to its vertex 


1 1 

V? 


= 4 Ax, 


1 u 
A’ 


(8.29) 

(8.30) 


(8.31) 


8.4.3 Algebraic spirals 

u = ^-^, \ = + (8.32) 

• m = 1: Archimedes Spiral. 

• m = 1: Hyperbolic or Reciprocal Spiral. 

• m = 2: Galilei’s Spiral. 

• m = Cotes’s Lituus. 

• m = ^: Fermat’s Spiral. 

8.4.4 Logarithmic or equiangular spiral 

• u = p = rsina, where a is the angle beween the curve and the radius vector. 


8.4.5 Sinusoidal spirals 

u = {cos k(j))^ , + (1 “ k^m?)u^. (8.33) 

pz 

If m = — 1 we obtain the Rhodenea u = — K-r. 

COS k(p 

k = 3 gives the Trifolium. 
k = 2 the Quadrifolium and 
k = ^ is the pedal of the Cardioid. 

If m = 1 we obtain the Epi-spirals u = cos kcj). 

If m = —^ we have u = {cos k(l))~k and ^ 

Some particlular sinusoidal spirals are: 

• m = 1, A: = ^. The Trisectrix of De Longe u = cos 2 ^ ^ j + f 

• m = 1, k = ^. Maclaurin’s Trisectrix u = cos f, ^ = 5 + |u^. 
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• m = 3,k = 1. The Cubic Duplicatrix u = cos^ (j), ^ = 9u3 — 8u^. 

• m = 2,k = 1. The Kampyle u = cos^ cj), ^ = Au — 3u^. 

• m = |, /c = 1. The Witch of Agnesi u = cos 2 (j)^ p = 4u3 — . 

• m = 2,k = 1. The Cruciform u = cos 20, ^ = 4 — 3u^. 

• m = —2, k = 1. The Oeuf Double u = cos“^ cj), ^ = 4n^ — 3u^. 

• m = —3, k = ^. Cayley’s Sextic u = —^ = us. 

8.4.6 Black holes 

As central orbits, the first two examples below are related by Arnold-Bohlin duality |39l 140] 
while the third is self-dual. They arise in the theory of null geodesics of spherically symmetric 
Ricci flat black hole metrics |41j . 

— coS ^-1 or u = ^+1 ^0 have ^ = 5 + These are particular null geodesics of a 

four-dimensional black hole. 

~ coS 2^-2 or u = we have ^ = | -|- These are particular null geodesics of a 

seven-dimensional black hole. 

If u = tanh or u = coth we have ^ = 2(1 -|-m^). These are particular null geodesics 

of a five-dimensional black hole. 

8.4.7 Other examples 

Freeth’s Nephroid is 

r= ^ 1 -F 2 sin^^ , 

whence 

1 ^ 9n 

~9 ~ T rt ) 

p 2 4 V 

8.4.8 Some hyperbolic spirals 

u = (tanh A;0)™', u = (coth /c0)"*, — = ™ ^ ^ (1 — 2m^k^)u^. (8.36) 

u = (cosh ^0)™", ^ =—m‘^k‘^u^~^ + (1 + w?‘k^)u^^ (8.37) 

pA 

n = (sinh A:0)”*, \ = m‘^k‘^u‘^~^ + {1 + m'^k‘^)u‘^. (8.38) 

pZ 

If m = 1 we get the two types of Poinsot Spirals for which 

= (A 2 + i) 2 ^^ 2 ^ (8.39) 

pZ 

respectivley. 


(8.34) 

(8.35) 
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8.4.9 Further examples 

• The astroid is 

1 . - 

— 1 

• The nehproid is 

J_ - 

4 — 16^2 

• The deltoid is 

1 _ 8ii2 

9v? — 1 

u = (tan/cc/))™", ^ _l_ {2m^k‘^ + l)u^. 

• Kappa ii = cot 4>, m = —1, /c = 1, ^ = 1 + 3u^ + 

• Moulin a vent u = cot 2(/), m = l,/c = 2, ^ = 4 + 9u^ + 4u^. 


(8.40) 

(8.41) 

(8.42) 

(8.43) 
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